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Abstract 

The genuine physical reasons explaining the delocalization effect of the Hubbard repulsion U are 
analyzed. First it is shown that always when this effect is observed, U acts on the background of a 
macroscopic degeneracy present in a multiband type of system. After this step I demonstrate that 
acting in such conditions, by strongly diminishing the double occupancy, U spreads out the con- 
tributions in the ground state wave function, hence strongly increases the one-particle localization 
length, consequently extends the one-particle behavior producing conditions for a delocalization 
effect. To be valuable, the reported results are presented in exact terms, being based on the first 
exact ground states deduced at half filling in two dimensions for a prototype two band system, the 
generic periodic Anderson model at finite value of the interaction. 
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I. INTRODUCTION 



The ability of the Hubbard repulsion U to drive transitions to insulating behavior is widely 
studied and knowi>i. However, in the recent years considerable interest is devoted to the 
study of the converse phenomenon, namely the delocalization effect produced by U—'^i^. This 
effect attracts much attention given by the unexpected nature of the process, since it can 
finally lead to an interaction driven metallic behavior starting from an insulator, instead 
to drive transitions to insulating behavior from a metal, as intuitively is expected. The 
interest for understanding the delocalization effect of U has been first set up by the observed 
conducting behavior in two dimensional (2D) disordered systems^ where the electrons are 
localized in the starting non-interacting phase by disorder. But, opposed to the random 
case, periodic external potential is also able to cause insulating behavior, hence can be as 
well used as a starting point in investigating the delocalization effect of U . On this line the 
ionic Hubbard model is investigated at half filling in different approximations, or numerical 
methods, starting from dynamical mean field theory^, it's cluster version^ or quantum Monte 
Carlo technique^, obtaining insulator to metal transition in 2D, even if several aspects of 
the results are under debate^. 

On the other side, the clear physical reason for the U produced delocalization is far to 
be properly understood. Up to this moment only correlation screening of the one-body 
potential responsible for the insulating phased has been invoked as an intuitive possibility 
for explanation. 

In order to understand the physical reasons for a possible insulator to metal transition 
caused by the Hubbard U we have to pay attention to two aspects which must concur 
in producing a conducting behavior starting from an insulator-: i) the one-particle states 
must become extended, and ii) the gap in the charge excitation spectrum must disappear. 
The correlation screening of the one-body potential can be in principle a physical insight 
in explaining point ii), but it does not explain how, given by f/, point i) becomes also 
satisfied, e.g. how the state becomes extended given by the Hubbard interaction in a genuine 
multielectronic case. The deep physical reasons of this last process are not clearly known, 
and in my knowledge have not been directly investigated so far at a genuine multielectronic 
level in conditions of interest here, e.g. around half filling. 

The purpose of the presented paper is to provide valuable information regarding this 
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issue. During the paper we aim to answer the question: Is there a direct evidence that t/, 
in the studied circumstances is able to extend one-particle states ? To obtain a valuable 
answer, the description is made here at exact level in 2D. Our strategy is to deduce in a non- 
approximated manner the exact ground states in conditions of interest. A such obtained state 
will be a genuine multielectronic state, hence usually, is not possible to discuss rigorously 
about one-particle states in its frame. But, based on it, one-particle behavior can be tested 
by calculating the expectation value of the long-range hopping type of correlations 

= ^^E^iW, r^r) = {l/2)Y,\{KMr,. + H.c)l (1) 

^ i=l (T 

where N\ represents the number of lattice sites, and creates an electron at site i with 
spin a. The study will be given in this paper for the non-disordered case in homogeneous 
conditions, hence r(r) = ri(r) holds. During such a study, the one-particle localization 
length £ can be also introduced by r(r) ~ exp(— |r|/£). In order to answer the question 
aimed to be investigated here, one simply must analyzes whether, or not, U is able to 
increase r(r) and I in the conditions of interest. 

In order to fix the frame of the study, we must clarify the nature of the circumstances of 
interest. On this line one observes that in both cases mentioned in the introductory part, 
the U driven delocalization effect emerges on the background of a macroscopic degeneracy. 
Indeed in the disordered cases the ground states are macroscopically degenerate^, while in 
the ionic Hubbard model case the effects are observed in the neighborhood of the point 
where the single occupancy at one-one site of the sublattices A and B have the same energy 
as an empty site in A and a double occupied site in B- . Furthermore, one knows that the 
ionic Hubbard model is of two band type^, while real materials, also disordered once, have 
a multiband structure which one usually projects in a few band picture^. This projection is 
stopped here for convenience at the two band level. 

Given by these observations, one realizes that in order to find the physical insights of 
interest here, one must considers a prototype two band model in 2D. The non-interacting 
case must be set in the position of a macroscopic degeneracy from the point of view of the 
one particle states such to find an insulator at f/ = 0. This can be achieved by choosing 
the Hamiltonian parameters such to provide a bare flat (e.g. macroscopically degenerate) 
lowest band, separated by a finite gap from the upper dispersive bare band, and fixing the 
filling of the system to half (e.g. lowest band: completely filled, upper band separated by a 
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gap: empty). After this step one turns on U, and analyzes its effect. 

This experience is quite interesting also because one knows that turning on the interaction 
in many-body systems with macroscopic degeneracy, one always finds a fascinating physics 
as for example in the case of the Laughlin state in two dimensions''^, Hubbard chains with 
bare degenerate bands'^, or geometrically frustrated spin systems^^. 

The prototype two band model chosen for the study in this paper will be the periodic 
Anderson model (PAM), simply because several technical aspects regarding its ground state 
construction are known^^ii^. The PAM has been constructed to be the simplest microscopic 
model used to investigate the characteristic properties of heavy fermion and intermediate- 
valence compounds containing elements with incompletely filled f-shells^*^, but serves as 
well, as a test prototype two band model used in various circumstances of interest^ii^ii^. 
The PAM consists of strongly correlated, almost localized f-electrons which experience an 
on-site Coulomb repulsion U > 0, and hybridize with a band of non-interacting conduction 
(d-) electrons. In contrast to the one dimensional Hubbard model, PAM is not integrable 
even in dimension D = 1. Therefore our knowledge of the physics described by the PAM 
at finite and nonzero U is almost exclusively obtained by approximate methods. Notable 
exceptions are the constructed exact ground states in dimensions D = 1,2, 3, at and below 
1/4, or at and above 3/4 fillingi^^^^iSO, e.g. n = N/{4:Na) e (0, 1/4] U [3/4, 1), where N 
is the number of electrons. The physically quite interesting region situated at and in the 
vicinity of half filling, strongly debated at the level of different approximations or numerical 
methoda^'^'^, lacks exact results for < t/ < oo. 

The used technique casts the Hamiltonian in a positive semidefinite form^^, and the fact 
that a such expression has the possible minimum eigenvalue zero, represents the route in 
obtaining the exact ground states^-. The method can be applied even in unexpected situ- 
ations in the context of exact solutions, as'^ D = 3, disordered and interacting systems in 
2D^, non-integrable Hubbard chains in external fields^, or stripe, checkerboard and droplet 
ground states'^ in 2D. However, given by strong correlation effects not treatable mathemat- 
ically up to this moment in the studied concentration region, it has not yet been applied for 
the PAM at finite U ^0 around half filling. 

The deduced results provide the first exact ground states for the 2D generic PAM at half 
filling and finite value of the interaction. This is achieved by the use of extended operators 
in constructing the ground state, technique which is described also for the first time in the 
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context of non-integrable systems in this paper. 

The comparison of ground states at f/ = and f/ > on the hne of ([1]) shows that the 
Hubbard repulsion, by strongly diminishing the double occupancy, spread out the contri- 
butions in the ground state, hence is able to increase by several order of magnitudes the 
one-particle localization length. Hence the one-particle behavior in the analyzed conditions 
is extended by U, which is a clear physical source of a delocalization effect. Consequently, 
the answer to the question we are looking for in this paper is: yes. I also note that one 
already knows that even away from half-filling, U is able to create in relatively similar con- 
ditions subtle correlated conducting phases^ii^*^. It is also known that insulator to metal 
transition can be driven as well by disorder—, the process emerging at half filling again from 
a strongly degenerate insulating phase. 

The remaining part of the paper is structured as follows. Section II. describes the model 
and its solution in the non-interacting case. Section HI. prepares the deduction by transform- 
ing the starting Hamiltonian in a positive semidefinite form. Section IV. describes the ex- 
tended operators that allow the construction of the exact ground states. Section V. presents 
the deduced exact ground state, characterizing as well the physical properties of the system, 
while Section VI. containing the summary and conclusions closes the presentation. 

II. THE MODEL USED 
A. The Hamiltonian 

I analyze a generic PAM Hamiltonian H = Hq + UUf defined on an A^a = L x L, 2D 
Bravais lattice with primitive vectors (x, y). The Hubbard term is Uf = n(-^h( ^, arbitrary 
U > is considered, and 

^0 = E { (E trA,A+r,,<r + Vodlk. + (2) 
i,o- ra 

E E Vr^'KA+r,. + H.C.) + EfflJ,^], 
r=x,yb,b'=dj;b^b' 

where in order to describe real systems, the hopping matrix element tj.^ covers by the index 
Yd nearest and next nearest neig hbors, Vo and , b,b' = dj, give the strengths of the 
on-site and nearest neighbor hybridizations, and Ef represents the on-site /-electron energy. 
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During calculations periodic boundary conditions are used, and the filling is fixed to half 
(e.g. n = 1/2, N = 2iVA). 



B. The [/ = non-interacting case 

In order to have a macroscopic degeneracy at one-particle level, one considers first the 
non-interacting f/ = case, and selects the Hq parameters such to obtain a lowest fiat band 
in the non-interacting band structure. For simplicity one considers the (x, y) symmetric case 
ti = tx = tyi ^2/2 = t2x = ^2y = iy±x/2, Vi = V^'^' = Vy'''' . To have a lowest fiat band in 
the non-interacting band structure one takes 

V,h = Voh, Ef = 2^-^-2t2, h>0. (3) 

Indeed, for conditions ([3]), Hq from ([2]) transformed in k space, can be diagonalized into the 
form 

cr k 

where Cl^^ represents the canonical Fermi operator of the z/th diagonalized band (here 
z/ = 1,2) of Ho- In the case of the two bands Ei,k = -{t2/2)[{V^ /V^) + 4] < E2,], = 
^i,k + 2t2[(^iA2)^ + (Ck + Vo/{2Vi)f], Ck = cos(xk) + cos(yk), are clearly separated, the 
lowest band is fiat, the band gap being A = 2t2[{Vo/ti)^ + {\Vo/2Vi\—2)^] ^ 0. Consequently, 
in the non-interacting case, given by the presence of the lowest bare fiat band, a macroscopic 
degeneracy is present, whose degree at the one particle level is 2N/<^. 
The canonical Fermi operators present in (jlj) are given by 

C'ik,. = ^[Vl-^k + stgnia^)^l + F^ flj, (5) 

where = Vo + 2V,C^, (3^ = V^Vta - t?/(4t2) - ti^ - t2Cl, and Fk = [1 + (ak//?k)T°-'- 
The ground state at f/ = and half filling (lowest band completely filled, upper band 
completely empty) becomes 

Na 

i^s)=nn^k.io), (6) 

o- k=l 



which represents a paramagnetic band insulator. The ground state energy corresponding to 

Denoting by (...)o ground state expectation values calculated in terms of (jS]), for each 
site i of the system one finds d{ = {n{i^n(^)Q ^ 0. Concerning r(r) from ([T]), using ([HD, 
b = d, f, one finds r(r) ~ exp(— |r|/£), where the associated one particle localization length 
i is of order unity in lattice constant units (for example, for r || x, one obtains ^/|x| = 0.5 
for the Hamiltonian parameters fixed in Fig.8.). Based on this result one concludes that the 
non-interacting system in the studied conditions represents a localized band insulator. 

III. THE HAMILTONIAN TRANSFORMED IN A POSITIVE SEMIDEFINITE 
FORM 

A. The transformed Hamiltonian 

In order to deduce the exact ground state in the interacting case, one first transforms the 
Hamiltonian H = Hq + Ulff, U > 0, given in (!2][3!) in a positive semi definite form. I note 
that this exact transformation is independent, and provides 

H = Ho + UUf = G + UUf + Eg{N), 

^ = E E E ^vA^^e, E,iN) = -K,N. (7) 

v=l,2 e=±l iGSi, 

To clarify Eq.([7]), the following explications must be given: In the first line of ([7]), G and 
UUf are positive semidefinite operators, while Eg{N) represents the ground state energy 
(see Sec.IV.), a number which has a given dependence on the Hamiltonian parameters, and 
whose expression will be clarified below. In order to obtain ([7]), first we divide the lattice 
in two sublattices. The sublattices are denoted by v = 1,2, and all lattice sites from the 
sublattice v are denoted by Sy. After this step one introduces the cell operators Ai ,, g) which 
are given by 
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Ai,D,e = Cl,/-Pl,A„,e + an,d-Di+r„,A„,e- (8) 

n=l 

Even if labelled by the site i, the Ai^^^e operators are not local, and collect contributions in 
the form of a linear combination, from five sites situated around the site i. These five sites 
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denoted by i + r„, n = 1,2,3,4,5, build up together the cell defined at the site i, where 
taken for increasing n, r„ has the values 0, — y,x, y, — x (see Fig.l). 

n=4 




i-y 

n=2 



FIG. 1: The cell defined at site i (thick line) containing 5 sites of the lattice (i, i it x, i it y) on 
which the operator A^^y^^ is defined. The sites of the sublattice v, i £ S^, are denoted by black 
dots, (x, y) are the primitive vectors of the lattice, and n = 1,2,. ..,5 represents the notation of 
sites inside the cell. 

After this step one presents below separately in details the numerical prefactors, and the 
operators present in 

B. The numerical prefactors of e 

The numerical coefficients ai j, an,d in dHD are related to the Hamiltonian parameters of 
-f^O) being given by the following system of equations 

Vf = al^aij, Vj'^ = alja34, V^'^ = aljai4, 

5 

Ef = laijl"^ - Kd, Kd = ^\an4\^. 

n=l 

It is important to note that the exact transformation presented in ([7j) works only if the 
,, matching conditions" ([9]) are satisfied. I also mention that the system of equations IQ, 
through the expression of Kd, fixes as well the value of the ground state energy Eg{N) in 
([7]). One obtains the system of equations ([9]) by 1) calculating Al^^Ai^y^^ based on flHfTTl) . 
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2) deducing G + Eg{N) in ([7]) starting from the previous result, and 3) finally equating to 
Hq from ([2]) the obtained expression. I note that the presence of /j^o- operator only in the 
middle of the cell operator Ai t, e in ([7]) (see also Fig.l), leads to the absence of the direct 
/-electron hopping in the transformed Hamiltonian. 

The matching conditions ([2]), in the (x, symmetric case specified above Eq.([3]), give 
solution only when ^ is satisfied. The solution is characterized by an>2,d = (^24, and one 
obtains 

ai,/-^i^e , «M-^y|e , - , (10) 

where is an arbitrary phase. 



C. The operators contained in A\^v,e 

The operators £)j_A„,e, -Pj,A„,e, e = ±1, are new genuine canonical Fermi operators substi- 
tuting the starting canonical Fermi operators Jj^o-, /j^o-, cr =t, i- Their expression is given by 
the following relations 

where for h = d (/), one has B = D (F), and A^, is an arbitrary complex number. 

Looking for the physical meaning of B operators, taking for example the b = d case, 
using the fermionic spin representation for 1/2 spin operators S^^ = (l/2)(nf^ —nf^), S^^ = 
dlAh S^,i = dlkh Sti = SI, ± tSl,, SI, = Sl,Sl, + il/2)iSX,S^, + ), and intro- 

ducing the states |A,A,e) = A^Ael^)' (-Di,A,e| A,A,e) = 1, one finds for arbitrary A 

^d,ilA,A,e) = ^|A,A,e), (nf^ ^ + ) I A,A,e) = |A,A,e)- (12) 

Since similar relations can also be written for the b = f case, one concludes that -Bj^;^^ 
introduces one spin 1/2 fermion {b = d, or b = f for B = D or B = F) on the site 
i. One further observes that at the level of expectation values, the spin direction and 
orientation are provided by A, and e. Indeed, for all a = x,y,z one has (Di,A,-i|5'^il-Di^A,-i) = 
-(A,A,il'5^,ilA,A,i), and 

(A,A,i|5,"i| A,A,i) = Y^^[^-."^e(^) + ^y,c.Imi>^) + sj-^^^], (13) 
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such that ((-Di,A,e|'S'd,i|-Di,A,e))^ = 1/4. These results show that in fact A fixes an arbitrary 
axis along which the spin is placed, while e = ±1 gives the spin direction along the fixed axis. 
Indeed, for example at A = 1 one has S'^i\Di x=i,e) = (e/2)|Di A=i,e), for A = i one obtains 
S'^i|A,A=j,e) = (e/2)|A,A=i,e), whilc A ^ oo corresponds to the SJilA.A.e) = (e/2)|A,A,e) 
case (e.g. for A = 1, A = z, and A oo, the spin axis becomes in order the x,y, and z axis). 

I further note that the relation ^ i^l^ _^ = dl ^dl^ holds for all A. It is also important 
to mention that in Eq.Q, after the transformation of the starting Hamiltonian defined by 
([2]), for f = 1, 2, one obtains arbitrary, but fixed Ai and A2. 

IV. DESCRIPTION OF THE MATHEMATICAL SOLUTION THAT LEADS TO 
THE GROUND STATE WAVE FUNCTION AT HALF FILLING AND U > 0. 

In ([7]) one has the Hamiltonian transformed in a positive semidefinite form which now 
must be analysed in order to deduce its ground state at half filling in the interacting U > 
case. As mentioned above ([2]), H has been defined on an L x L Bravais lattice, and obtaining 
the transformed form in ([7]), the lattice has been divided in two sublattices f = 1, 2 containing 
each L^/2 sites. 

In deducing the ground state at particles we have to find a wave vector \(j){N)) = K'^lO), 
where introduces electrons into the system, and based on (I7l) 

G\^(N))=Uf\^iN)) = (14) 

holds. In such conditions \(j){N)) represents the ground state of the A^-electron system 
described by in (^^, the ground state energy being Eg{N). In order to prove f |T4l) . one 
shows {Ai^y^e, = for the first part, and demonstrate that \(f){N)) not contains double 
/-occupancy for the second. 

Below I describe the K'' operator which will provide the ground state at = 2L^ (half 
filling). Since in this concentration region strong correlation effects are present, K'^ has a 
relatively complicated form. Because of this fact, in describing first I present it, describe 
it in details, give its expression, and explain how it was deduced. 
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A. The mathematical key of the solution 



Taking |0) as the bare vacuum, one considers the state 

2 

l*e^,i^) = n n (^i,i..=l,A„,e>i,.=4.A.,eJ|0), (15) 

v=i ieSv 

where W^^^ is the broken hne operator defined for the sublattice v, labelled by the site i 
of Sy, breaking direction at the site i Ti = r, spin direction given by A„, and spin orientation 
specified by Cy. 

One arrives to the broken line operator notion by looking for the desired 2L^ number of 
linearly independent operators suitable in constructing the ground state at n = 1/2. If we 
would like to use an operator of the type Cl^.^ (or its Fourier transform in r space) from 
f lS|^ for this job, we shortly realize that only such operators are available in constructing 
the ground state aX U > 0, namely those with fixed a, since mixing the spin index, auto- 
matically double occupancy occurs which raises the energy. Hence, the expression of the 
interacting ground state requires new operators W"^ which must be find. The search for these 
leads first to operators constructed along lines in y + x direction. These resemble to the 
operators presented in Fig. 7, but their number is of order L, hence is small. At this point 
one remembers that by breaking a straight line, a site (e.g. the breaking point) becomes 
attached to the line, the number of sites being L^. Hence in this manner is possible to reach 
the number of desired operators. The presence of the cusp point on the line complicates a 
little the expression of the operators, but this is in fact the way one reaches the broken line 
operator notion. 

Since the structure of the operator is not transparent in (fT5|) . I explain it in details 
below. 

B. Straight lines and broken lines 

Connecting the sites of a given sublattice in y + x direction, n e [1, L/2] different straight 
lines /y+x,t;,n can be drawn, each containing L sites, see Fig.2. For simplicity one considers 
below even L. 

Now one breaks the straight lines at an arbitrary site i G ly+x,v,n- In order to be this 
possible, one first defines 4 breaking directions denoted by the index r as shown in Fig.3. 
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FIG. 2: An 8x8 lattice (e.g. L = 8 holds), containing one straight line ^y+x,'y=i,n) (thick line). 
The direction of the line is x + y, it is contained in the sublattice v = 1 denoted by full circles, 
the line contains L = 8 sites, and note that periodic boundary conditions are considered. (x,y) are 
the primitive vectors of the 2D Bravais lattice. In the x + y direction L/2 = 4 such lines can be 
plotted. 

If one fixes the breaking direction at i by ti = 1, one finds the ,, broken line" as presented 

V > ^ < 

»2 « 

FIG. 3: The four possible breaking directions denoted by the index r. 

in Fig.4. This is simply obtained by taking the straight line /y+x,'y=i,n crossing the site i, 
and breaking the line at the site i in the breaking direction Ti — 1. Given by the periodic 
boundary conditions, a broken line besides its first cusp at i, and — 1 (obtained by breaking 
the starting straight line at i), will have a second cusp point at j = i+(L/2)(x+y), and tj = 3 
(see Fig.4), which is unambiguously determined by (i,Ti). Hence the obtained broken lines 
can be denoted by lv{i,Ti — 1), and their number for fixed v becomes Consequently, 
the broken line /^(i, Ti = 1) is a continuous line in the sublattice v, containing L sites from 
Sy and two cusps as presented in Fig.4, the cusp characterized by the breaking direction 
T — 1 being present at the site i. 
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FIG. 4: A broken line lv{hTi = 1) (thick line) defined at the site i in a 8 x 8 lattice. Ti = 1) 
connects sites of the v sublattice (full dots). Note that such broken lines can be plotted for 
a fixed V. 

C. The operators 

The l^Jii-A„e„ operator represents a linear combination of operators acting on i' sites 
placed in an extended region surrounding the broken line Ti = r) as presented in Figs. 5- 
7 for r = 1. In these figures, the black, open, and gray circles plotted at a site i' are 

representing in order A',A,,e, = ^^f' A^A„,e„ + 4^',A«,e„' ^f',A„,e„. and ^!,A„e„ operators. 

The coefficients present in Figs. 5-7 are numerical prefactors. In order to obtain the 
mathematical expression of the l^JiT-=iA^e„ operator, one simply adds all contributions 
from Figs. 5-7 as follows 

^v,i,T=l,Xv,ev ~ ^-^-2y,A^,e^ + A^_y,A„,e„ + ^^i+y,K,ev + /^-^+2y,A„,e^ 

M 



^ Z^V-^n-^i+n(x+y),A„,e„ ^ -^n^ i+n(x+y),A„,e„^ 
n=0 

M-1 

+ (^n-^i+n(y-x),A„,e„ + ^n-^+n(x+y),A„,e„) 

n=l 

M-1 

+ ^2 ^"-^i+nx+(n-l)y,A^,e^ + X/ ^"-^i-nx+(n-l)y,A„ 
n=l n=l 
M M 

+ ^?/n-Di+(n-l)x+ny,A„,e„ + X] A 



ra=2 



)x+ny,A„,e„ ' / ^ i-(n-l)x+ny,A„,e„ 
n=2 



' J-y,A^,e„ j+y,A„,e„ 



^ " j-2y,A„,e„ 



+ «l^+y+x,A„,e„ + «l^Vy-x,A„,e„ + «oi^V2y,A„,e„ > 
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(16) 



where M = L/2, j = i + (L/2)(x + y), and effectuating the sums over the index n in (fT6l) . 
the periodic boundary conditions must be taken into account. 




FIG. 5: The first cusp region of W^^-^ j ^^-^ constructed along the broken hne (thick line) 
labelled by (i, Tj = 1). White, gray and full circles represent in order l)^, F^, and -P^ operators, 
the coefficients {xi,yi,e, etc.) are the numerical prefactors. 
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FIG. 6: The second cusp region of 

=1 i T=i A„ El,- '^^^ thick line is a part of the broken line. For 

notations see Fig. 5. 



A similar procedure based on Ti = 4, tj = 2 leads to the broken line used in constructing 
W}- ^. ^ ^ , which can be obtained from Wl- , ^ „ by the interchange x ^ y. 
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FIG. 7: Inter-cusps (straight) region of W^J^j^ ; ^^^^ ^ e • "^^^ thick hne is a part of the broken hne. 
Notations as in Fig. 5. 

D. The numerical prefactors of the operators 



The numerical prefactors present in Figs. 5-7 and Eqs.( |T6|) are given as follows 
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(17) 



aij aij 

where yo is arbitrary, and 

The deduction procedure of (IT7|) is the following one: 1) One fixes a lyj^; operator 
labelled by the fixed site i of the sublattice v as shown by Figs. 5-7. 2) After this step one 
takes all operators Ai/^^/^e' of the Hamiltonian from Eq.([7]) which share common sites with 
^viiXyCy^ 3) calculate the anti-commutators { Aiy^e' , ^vHXyCy } ' 4) require zero 
value for these, e.g. 

{ A'y,e' , W^i,M,A„,e„ } = 0, (18) 

for all possible values of all indices. These conditions provide a linear system of equations, 
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whose linearly independent terms build up the following system of equations: 
In the region of the lower cusp one has 

OljC + «2,d|/0 = 0, 

aijii + a2,de + 2yoa24 = 0, 
a2,d{yo + e) + ai,dxi + aijx( = 0, 
a2,dXo + ai^dVo = 0, 
ai,dxi + aijx{ = 0, 

ai,d^i + (^ij^i = 0> 
eoi.d + a2,d{x^ + 2xf) = 0, 
a2,d{xo + xf) + ai^yi = 0, yi = -yo, 
a2,d{^i + ^2) + «i,d2/2 = 0, y2 = yo: 
a2,d{^2 + xf) + ai^dyz = 0, ya = -yo, 
a2,d{xi + 3:4) + ai^y^ = 0, ^4 = Z/o- 
For the inter-cusp region presented in Fig. 7 one obtains for 4 < m < (M — 

+ aijxl^ = 0, 
a2,d(a:^ + = -ai,d2/m+i, 

= {-iTyo. 

The upper cusp region provides the equations 

aij5 + a2,dl ^ -2a2,dZ/M-i, 

2a2,dyM + a2,d{P + 7) + ^i.^x^f + Qija^t = 0, 

aijtto + a2,d/3 = 0, 

a2,d(7 + Z/m) + aijXM^i = -2a2,dyM-i - ai^da^tf-i, 

a2,d{P + yM) = -aijai, 

a2,diyM + yu-i) = 0, 

«2,da^M + cn,d7 = ~2a2,da;^_i, 

a2,da:;M + = 0, 

0'2,dXM + 0>l,dyM — —0'2,dX'l^_^ 
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The system of equations Eqs. fll9ll20|[2Tl) must be solved considering yo a given, but arbitrary 
quantity. The solution is presented in Eq. (|T7|) . 

E. The properties of the operators 

The properties related to the constructed i ^ a„ e„ operators are the following ones: 

a) Since different possible sites i € S*^ are present in a given sublattice v, the operator 
4>e:iv) = niG5.(^i,r=i,A„,e„W^Ii,r=4,A„,eJ mtroduces elcctrous into the system. Hence, 
based on [nD=i 2 *^ej(^)]' ^^e vector from f[T^ becomes 

\'^eli') = [U^e:iv)m, (22) 

v=l,2 

and has 2L^ electrons, e.g. is defined at half filling. 

b) The operators i r a„ e„ Present in (!22l) are constructed from 2L^ canonical Fermi 
operators (-Di A„,e„5 A,A„,e„)^5 and are linearly independent, since each of them has at least 
one new segment^ in its construction. 

c) Since Wj;^^ g^PVj; = holds, the operators Wu,i,r,A„,e„ are fermionic opera- 
tors, but not satisfy canonical anti-commutation rules, e.g. {Wv^i^r,K,e^,Wj^, ^, ^ e' } ^ 

5v,-u'5i,i'^r,T"5A,A'5e,e'- 

d) Given by the construction of W^^^^^ operators [e.g. f|T8l) ]. one has 

G\^eli) = 0' (23) 

hence (1151) give the minimum energy eigenvalue (e.g. zero) of the first term of the studied 
Hamiltonian from (jT]). 

e) A quite important property of the W^J;^;^ g operators is that they introduce / elec- 
trons only in the sublattice v. This can be seen for example from Figs. (5-7) which best 
present the W"^ operators, by taking into account that F"^ operators (containing the 
operators) act only on sites depicted in these figures by gray and full (black) circles. Conse- 
quently, the 0ej(f) operator introduces / electrons only in the sublattice v, and more than 
this, because of -fl|A„,e„-^!A„,e„ — 0; vector 0g"(t>)|O) will contain only single / occupancy 
in the sublattice v, while the other sublattice v' ^ v will contain only d electrons. Conclud- 
ing this point, the contributions created by 0ej(f)|O) will have i) in 5^ empty sites, single d 



17 



or / occupancy, and double occupancy of the form dIx^ e^^iXy e„ sites i E S^, while ii) 
in the sublattice S'^,/, v' ^ v empty sites, and single d occupancy can appear. 

f) Because of the property mentioned at point e), l^^^l'^^) from (!22|) will contain F^Xi ei 
operators in the sublattice f = 1, i G Si; Pj^^ operators in the sublattice v = 2, ] E S2, 
hence double / occupancy is not present, consequently one has 

UUf\^'el'X')=0. (24) 

In conclusion, (fT5l) gives the minimum energy eigenvalue (e.g. zero) also for the second term 
of the studied Hamiltonian from ([7]). Note that even if Ai, A2 are arbitrary but fixed in ([7]) 
hence also in fl23|24p . the parameters 61,62 can be arbitrary chosen. 

g) As seen in Figs. 4-7 the operators are extended. The extended nature has several 
aspects that has to be separately stressed. 1) First, it is not possible to find 2L^ solutions 
of (fT8|) for W'^ operators such to obtain W'^ enclosed in a cluster smaller than the zise of 
the system, and to preserve the property ([23]). 2) By plotting all W"^ operators present for 
a fixed L one realizes that for two arbitrary sites (ji, j2) of the system there is at least one 
W'^ operator containing both sites, hence representing a direct path between ji and }2- 

V. THE GROUND STATE AT ?7 > 
A. The ground state at half filling 

The results presented in Sec. IV. E. [see ( 123|2^ ] show that the ground state at half filling 
of H is provided by the I^E'e^ '^^) vectors. Since as specified in Sec. IV. E, 61,62 = ±1 can be 
arbitrarily taken, the ground state at half filling becomes 

\^,{2L')) = \^'^''^)= J2 I'^el'i'), (25) 

ei,e2=±l 

the ground state energy being Eg(2L'^) = —2L'^Kd. Using fpiTUl) one finds Kd = 
it2/2)[{Vo/Vif + 4], hence Eg{2L'^) = ^°(2L2) holds for arbitrary U > 0. I stress here 
that since double /-occupancy is not present in the interacting ground state, only the pres- 
ence of f/ > is required for |\E'g(2L^)), hence f l25|) is the ground state for arbitrary U > 0. 

Defining 0'^" = Y^e^^ey^ ^^e vector ^'^''lO) represents a singlet state for arbitrary A^,. 
Since the ground state (125|) becomes now |\E'g(2L^)) = Jl^r |0)) represents a singlet, 
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non- magnetic state. Given by the arbitrary nature of the Ai, A2 coefficients, one could 
write the ground state fl25l) also in the form \'^g{2L^)) = XIai A2 ^^i.^al^^^''^^) where cai,a2 
are arbitrary numerical coefficients. But this expression not emphasize a supplementary 
degeneracy given by the parameters Ai, A2, only shows that the spin axis for both v = 1,2 
can be arbitrarily chosen. 




J t 



2 4 6 r 

FIG. 8: Ground state expectation value of f-hopping correlations r(r) at [/ = (dashed line) and 
[7/0 (continuous line) for a 12 x 12 system at ti/t2 = 10, ^2/^1 = 0.5, Vi/Vq = 0.1. r || (x + y) is 
in |x + y| units. Inset shows the log-log plot at U 7^ 0. For example, at r = 4, Tij/Tij=q = 2 • 10^ 
holds^l. The result is qualitatively the same for other directions and other Vi/Vq ratios as well. 

B. The one particle localization length in the interacting case 

The obtained result in (1251) shows that turning on an arbitrary U > 0, the ground state 
at n = 1/2 has the property d( = on each site i, where d( = {h(-^h(^), and (....) denotes 
ground state expectation value calculated in terms of |\&g(2L^)). Consequently, since d( > 
at f/ = 0, the Hubbard repulsion, by strongly diminishing the f-electron double occupancy, 
spread out the contributions in the ground state wave function (e.g. the ground state in exact 
terms can be now given only in terms of operators which extend along the whole system). 
The aspects of this fact have been underlined in the point g) of Sec.IV.E, these properties 
are missing in the non-interacting case, and are introduced by the Hubbard U > 0. As 
a consequence, if one fixes two arbitrary sites (jnjjm) of the system, there is at least one 
^i!irA^e„ Operator which contains both sites, hence represents a direct path between j„, 
and j„, e.g. gives {Wv^,T^,,eJ^^^„kn^,^\Wv^,r^,,eJ / {W^^,T^,,eJW^ ~ 0{1/L), where 

\W^,i,T,x,,eJ = W^J,i,r,A„,ejO)' ^nd r = |j„ - j^l is arbitrary28. 
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This property produces a clear increase in the one particle localization length i. I show 
this aspect by presenting an exemplification for L = 12 in Fig. 8, where r(r) from ([1]) 
calculated in terms of |\l/g(2L^)) is shown. This emphasizes that r(r) has no more strong 
exponential decay as in the f/ = non-interacting case, and can increase even 5 order of 
magnitudes given by the presence of a non-zero U . The i value corresponding to r(r) for 
the interacting case in Fig. 8 (see inset) exceeds the sample size, while at [/ = had the 
value £/|x + y| ~ 0.5. The delocalization effect produced by the Hubbard f/ > is clearly 
seen. 

The reason for this behavior is that U throw the electrons about in order to diminish 
considerably the double occupancy. This is why W'^ becomes extended along the whole 
system. In other words, it is not possible to construct at half filling in the presence of 
f/ > W'^ operators with the properties presented in Sec.IV.E., enclosed in a cluster with 
fixed extension smaller than the sample size. Concerning the needed degeneracy, I note the 
following: In order to construct the resulting extended W'^ operator from ( fT6|) . the system 
must combine a huge number of linearly independent contributions. If all these contributions 
are belong to the same energy (e.g. the starting point is strongly degenerate), the creation 
of the interacting ground state containing the operators is without cost of energy, hence 
it is easier to be effectuated. This is the reason why macroscopic degeneracy is important 
for the process. 

C. The compensation effect in the PAM 

I add at the end an observation regarding the compensation effect in the PAM which is also 
visible from the obtained solution. In the PAM, the /-electron moments are compensated in 
the singlet ground state at n = 1/2. Starting from the paper of Nozieres^S^ the global or local 
nature of this compensation effect, especially in the large V limit, is strongly debated. The 
presented results show that for arbitrary large repulsive f/ values the compensation can be 
manifestly global, being given by a subtle superposition effect of arbitrary spin orientations, 
involving both b = d, f electrons. 
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VI. SUMMARY AND CONCLUSIONS 



In the study of the delocahzation effect of the Hubbard interaction U > 0, two aspects 
have to be analyzed in order to understand the physical reasons of this process: i) how U 
extends the one particle states, and ii) how, given by U, the gap in the charge excitation 
spectrum disappears. For the point ii) the results deduced for the ionic Hubbard model at 
half filling^i^ii^ suggest that the correlation screening of the one-body potential responsible for 
the insulating phase can be invoked as an intuitive possibility for explanation. This aspect 
however not emphasizes clearly the physical reasons for point i), namely, how the Hubbard 
repulsion extends the one-particle behavior in a genuine multielectronic and interacting case. 

The presented paper provide physical insights regarding this aspect i). First one shows 
that always, when the delocahzation effect of the Hubbard repulsion is observed, U acts 
on the background of a macroscopic degeneracy in a multiband system. After this step 
the delocahzation effect is analyzed produced on a band insulator starting point taken for 
the non-interacting case. For this one considers a prototype two band model at half filling, 
which at f/ = has a completely filled lowest band, separated by a gap from the empty 
upper band. In order to have the macroscopic degeneracy present, the lowest band has been 
chosen flat. After this step U has been turned on, and its effect analyzed. The extended 
nature of the one-particle behavior was tested by the ground state expectation value of the 
long-range hopping terms. 

This analyzis is made in the presented paper based on the periodic Anderson model 
(PAM) in two dimensions and half filling as a prototype model. In order to provide valuable 
results, the study is presented at exact level. The PAM has been chosen since several aspects 
are known regarding the construction of its ground states^^i^, but, exact ground states for 
it in 2D, interacting case, and half filling, are deduced here for the first time. This is 
achieved by the use of extended operators, the technique being described in the context of 
non-integrable models, again for the first time in this paper. 

The deduced results show that acting on the background of a macroscopic degeneracy, 
the Hubbard repulsion, by diminishing the double occupancy, spread out the contributions 
in the ground state wave function. By this, is able to increase considerably the one-particle 
localization length, hence the one-particle behavior becomes extended, which represents a 
clear physical source of the delocahzation effect. 
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The study of the point ii) mentioned previously, and related to the gap disappearance 
needs the deduction of the low laying excitation spectrum in the presented frame. The 
description of this case, at the level of the same quality, remains a challenging subject for 
future investigations. 
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